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Abstract: 

We obtain several higher order periodic solutions of a Coupled (f)^ model in terms of Lame polynomials 
of order one and two. These solutions are unusual in the sense that while they are the solutions of the 
coupled problem, they are not the solutions of the uncoupled problem. We also obtain exact solutions of 
coupled (p^-(/)'^ models, both when the 4>^ potential corresponds to a first order (asymmetric double well) 
or a second order (symmetric double well) transition. 
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1 Introduction 

Coupled triple well or (j)^ models arise in the context of many first order structural phase transitions. 
There exist analogous coupled models in field theoretical contexts [31 IH E] • Specifically, when a first order 
transition is driven by two primary order parameters, the free energy should be expanded to sixth order 
in both order parameters with a bi-quadratic (or possibly other symmetry allowed) coupling. In a recent 
publication we obtained a large number of periodic solutions of a coupled (j)^ model with bi-quadratic 
coupling [B]. All these solutions had the feature that in the uncoupled limit, they reduce to the well known 
solutions of the corresponding uncoupled cj)^ problem. The purpose of this paper is to point out that this 
coupled model has truly novel solutions in terms of Lame polynomials of order one and two [T', 'S], provided 
we add (symmetry allowed) quartic-quadratic and quadratic-quartic couplings to the coupled (j)^ model 
considered earlier [6] . It may be noted that these solutions exist only because of the coupling between the 
two fields (up to sixth order). In other words, while the Lame polynomials of order one and two are the 
solutions of the coupled problem, they are not the solutions of the uncoupled problem. 

In this context, we note that in recent publications we have obtained a large number of periodic 
solutions, in terms of Lame polynomials of order one and two, of coupled (j)'^ problems, both when the 
(j/^ potential corresponds to the second as well as to the first order transition [9l [TOl [TT] . The obvious 
question then is whether one can also obtain exact solutions of the coupled problems, both when 

the (/>^ potential corresponds to the second as well as to the first order transition. Interestingly, examples 
of these situations occur in condensed matter. The face-centered cubic to monoclinic transition in Pu 
involves an intermediate phase hexagonal to monoclinic transition with a free energy modeled by the 
coupled (/)®-asymmetric model [12]. Similarly, the coexistence of face-centered cubic, body-centered 
cubic and hexagonal close-packed structures in cobalt [13] as well as Fe and Tl |14] is modeled by the same 
free energy. An example of coupled (/)^-symmetric 0^ model with a bi-quadratic coupling is the triggered 
ferroelectric transition [151 116j. Thus, another important purpose of this paper is to obtain exact solutions 
of the coupled cfP-cf)'^ models with bi-quadratic coupling, both when the 0^ potential corresponds to the 
first as well as to the second order transition. 

The plan of the paper is the following. In Sec. II we provide novel periodic as well as the corresponding 
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hyperbolic solutions in terms of Lame polynomials of order one for the coupled (p model with an explicit 
bi-quadratic as well as quadratic-quartic and quartic-quadratic couplings. In Sec. Ill we provide novel 
periodic as well as the corresponding hyperbolic solutions in terms of Lame polynomials of order two for 
the same coupled (p^ model. In Sec. IV we provide the exact solutions of the coupled cp^ — cp'^ problem with 
bi-quadratic coupling in case the cp"^ potential corresponds to either the first or the second order transition. 
Finally, in Sec. V we conclude with summary and possible extensions. 

2 The Coupled 0^ Model and Solutions in Terms of Lame Polynomials 
of Order One 

In P we had considered the following coupled (p^ model, with a bi-quadratic coupling, in one dimension 
with the potential 

Vi<P, ^) = [^<P' - \<P' + |/) + (f - + 1^^) + ^0^^ . (1) 
We now show that in case we add the following quartic-quadratic and quadratic-quartic coupling terms 

^' = |</'V + ^<^V, (2) 

to the potential ([T]), then in addition to the solutions obtained in [B], there exist truly novel solutions in 
terms of Lame polynomials of order one and two to this coupled problem. Here ai,2) &i,25 ci^2) d^ e and / 
are material (or system) dependentparameters; and 'ip are scalar fields. From stability considerations we 
shall always take ci,C2 > 0. Further, since we are interested in a model for first order transition, we shall 
take 6i, 62 > 0. As far as ai, 02 are concerned, their sign is arbitrary and the shape of the potential depends 
on the ratio 6f/4aiCi and 62/40202. In particular, in the decoupled limit (i.e. d = e = / = 0), it is easily 
shown that as long as 4aiCi > 6f , the potential has a minimum at i?i> = \yi\ I18j . In the case 4aiCi = 
apart from the minimum at = one now has points of infiection at = hi/ 2ci . As ai decreases further 
so that 4aiCi < h\ < (16/3)aiCi, one finds that while (/> = is still the absolute minimum, one now has 
two local minima and two maxima at 



At the special value of (16/3)oiCi = 6f one has three degenerate minima at cj) = and at (/)min given by 
Eq. (I2D; for relevant figures, see |17t I18j. This is the point of first order transition. As ai decreases further 
so that < (16/3)aiCi < 6f , then the roles of the minima are reversed, now i;^ = is the local minimum 
while (p'^in ^ given by Eq. ^ are the two degenerate absolute minima and cp'^^^ are the two maxima. 
Finally as oi < 0, the potential has two absolute minima at (pmin ^ given by Eq. while (j) = is now 
the sole maximum. This picture continues to persist, for arbitrarily large and negative ai. Throughout 
this paper, we shall term the point where bf = (16/3)aici as the point of first order transition, i.e. the 
point with T = Tc- On the other hand, the point where local structure (i.e. local minima) starts growing 
with decreasing temperature, i.e. h\ = 4aiCi, will be termed as T = Tp. Note that for T > Tp there is 
only a global minimum at = and there are no other extrema. Thus the region where 6f > (16/3)aiCi 
corresponds to T < Tc while the region with 4aici < 6f < (16/3)aici corresponds to Tc < T < Tp. A 
similar analysis is also true for the potential in ^p with oi, 6i, ci being replaced by 02, 62, C2, respectively. 
The (static) equations of motion which follow from Eqs. ([T]) and ^ are 

^ = a2^- b2^^ + C2^^ + di/jcj)^ + + /<^V^ • (4) 

These coupled equations have thirteen distinct periodic (elliptic function) solutions, i.e. five "bright- 
bright", three "bright-dark" and five "dark-dark" solutions which have already been discussed in [6] (in 
case e = / = 0). In particular, there are five solutions below the transition temperature Tc, four at Tc, one 
above Tc (i.e. Tc < T < Tp), and three in the mixed phase in the sense that while one of the field is above 
Tc, the other one is below Tc- The latter situation is akin to the one found in multiferroic materials where 
one transition (i.e. antiferromagnetic) takes place at a higher temperature than the other transition (e.g. 
ferroelectric) or vice versa p!9]. It is worth pointing out that in turn in the single soliton limit, these lead 
to eight distinct coupled (hyperbolic) soliton solutions. In particular, one obtains three solutions below the 
transition temperature Tc, two at Tc, one above Tc and two in the mixed phase in the sense that while one 
of the field is above Tc, the other one is below Tc. 

We now show that apart from these solutions, we also have rather unusual solutions in terms of Lame 
polynomials of order one (and two discussed in the next section) which we now discuss one by one. Since 



there are three Lame polynomials of order one (i.e. sn, cn, dn) and since the field equations are essentially 
symmetric in (p and ^/J, we expect six independent solutions to the coupled field equations in terms of Lame 
polynomials of order one. In particular, we first show that there are three periodic bright-bright, two 
periodic dark-bright and one periodic dark-dark soliton solutions in terms of Lame polynomials of order 
one, which in turn lead to one bright-bright, one dark-dark and one dark-bright hyperbolic soliton solution. 

2.1 Solution I 

We look for the most general solutions to the coupled Eqs. (jH) in terms of the Jacobi elliptic functions 
sn(x,m), cn(x,m) and dn(x,m) [7] where the modulus m = k"^. It is easily shown that 

(j) = Asn{Dx + XQ,m) , ip = Bsn{Dx + XQ,m) , (5) 

is an exact dark-dark periodic solution to the coupled Eqs. (jU provided the following six coupled equations 
are satisfied 

ai = -{l + m)D^ , (6) 

- biA^ + dB^ = 2mD^ , (7) 
4 + eA^B^ + ^ , (8) 

a2 = -(,l + m)D^ , (9) 

- 62^^ + dA'^ = 2mD'^ , (10) 

C2B^ + fA^B^ + ^A^ = 0. (11) 

Here, A and B denote the amplitudes of the kink lattice, D is an inverse characteristic length while xq 
is the (arbitrary) location of the kink. Three of these equations determine the three unknowns A,B,D 
while the other three equations give three constraints between the nine parameters ai^2i bi,2,ci^2Td, e, f. In 
particular, we find that the solution exists only if ai < 0, 02 < 0, e < 0, / < 0. We obtain 

n2_ ,2 _ 2mD^{d + b2) ^2 _ {d + bi)A^ 

"(1 + m)' ^ " d^-b,b2 ' " (d + 62) ' ^ ' 
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while the three constraints are 

ai = 02 < , (4ciC2 - \e\\f\){d + 61) = 2{e^ + 2|/|ci)(d + 62) , 

(4ciC2 - |e||/|)2 = 4(e2 + 2\f\ci)if + 2|e|c2) . (13) 

In the hmit of m = 1, the periodic solution ([5]) goes over to the hyperbolic dark-dark soliton solution 

(j) = A tanh(_Dx + xq) , = B tanh(I?x + xq) , (14) 

provided the constraints ()12p and (I13p with m = 1 are satisfied. 

2.2 Solution II 

It is easy to show that 

4) = Acii[Dx + XQ^m) , ip = Bcn{Dx + xo,m) , (15) 

is an exact bright-bright periodic solution to the coupled Eqs. (jH) provided six coupled equations similar 
to Eqs. dni) to (jlip are satisfied. Three of these equations again determine the three unknowns A,B,D 
while the other three equations give three constraints between the nine parameters ai^2, bi^2,ci^2,d, e, f. In 
particular, we obtain 

2 ai ,2 2m{d + b2)D'' „2 (^1 + d)A^ 

^ = 7^ TT ' «i = «2 , A = — — , B = , (16) 

(2m -1) (6162 -a^) [02 + a) 

while the remaining two constraints are again given by Eq. (|13p . Note that oi = 02 > (<)Oifm> (<)l/2. 
In the limit of m = 1, the periodic solution (jlSp goes over to the hyperbolic bright-bright solution 



(p = Asech^Dx + xq) , ip = Bsech.{Dx + xq) , (17) 
provided the constraints (fT3|) and ([16]) with m = 1 are satisfied. 



2.3 Solution III 

Yet another bright-bright periodic soliton solution is 

(p = Adn{Dx + XQ,m) , ip = Bdn{Dx + xo,m) , (18) 



provided six coupled equations similar to Eqs. ([6|) to pT|) are satisfied. Three of these equations deter- 
mine the three unknowns A,B,D while the other three equations give three constraints between the nine 
parameters ai^2,bi^2,ci^2,d,e, f. In particular, we obtain 

p2_ «i „ -„ ,2_ 2(rf + 62)ZJ2 p2_ ibl+d)A^ 

^ ~J2^)' ^ - {b^h2-<P) ' ^ - (62 + d) ' 

while the remaining two constraints are again given by Eq. (jl3p . 

In the limit of m = 1, the periodic solution (|18p again goes over to the hyperbolic bright-bright soliton 

solution p!7|) . 

Note that while for the solution d? > 61621 for the solutions and (fTHj) . its the other way around, 
i.e. d? < 6162. 

2.4 Solution IV 

Yet another bright-bright periodic soliton solution is 

(p = A^/mcn{Dx + XQ,m) , ^p = Bdn{Dx + XQ,m) , (20) 
provided the following six coupled equations are satisfied 

ai + ^(1 - m)'^B^ + d{l - m)B'^ = (2m - 1)D'^ , (21) 
61^2 _^ ^^2 _ (1 _ rn)eA^B'^ - (1 - m)fB'^ = , (22) 
cxA^ + €^2^2 ^ ^ Q ^ (23) 



02 + - "^)^^^ - (1 - ra)dA^ = (2 - m)!)^ ^ (24) 
62^^ + dJ^ + (1 - m)fA'^B^ + (1 - m)e^^ = 2D2 ^ (25) 
C2S^ + + ^ Q _ (26) 

Three of these equations determine the three unknowns A,B,D while the other three equations give three 
constraints between the nine parameters ai^2ybi^2yCi^2yd,e, f. In particular, this solution is also valid only 
if e < 0, / < 0. Further, two of the relations are given by 

(4ciC2 - |e||/|)2 = 4(e2 + 2|/|ci)(/2 + 2\e\c2) , (4ciC2 - \e\\f\)A^ = 2{f + 2\e\c2)B^ . (27) 



In the limit of m = 1, the periodic solution ()19p again goes over to the hyperbolic bright-bright soliton 
solution p!7|) . 

2.5 Solution V 

In addition there are two dark-bright periodic soliton solutions. One of them is 

(j) = Asn{Dx + xo,m) , tp = Bcn{Dx + xo,m) , (28) 
provided the following six coupled equations are satisfied 

ai + ^B^ + dB'^ = -{1 + m)D^ , (29) 

- biA^ - dB^ + eA^B"^ - fB^ = 2mD^ , (30) 

ciA^ + ^B^ = eA^B^ , (31) 

02 + + dA^ = (2m - 1)1)2 ^ (32) 

62-6^ + dA^ - fA^B^ + eA^ = 2mD^ , (33) 

C2B^ + = fA^B^ . (34) 

Three of these equations determine the three unknowns A,B,D while the other three equations give three 
constraints between the nine parameters ai,2i ^i,2) ci^2) e, /. For example, it is easily shown that 



^ 61 ± J62 - 4aici - 4(1 - m)Z)2ci ^ b2 ± Jbl - Aa2C2 - 4:D^C2 

s;; ■ ^ - 2c, ■ '='^> 

Further, unlike the previous four solutions, this solution exists only if e > 0, / > and two of the constraints 
are given by 

(4ciC2 - ef)A^ = 2(2ec2 - f)B^ , (4ciC2 - eff = 4(2ec2 - f){2fci - e^) , (36) 

while the other two constraints are 

2 _ 61^2 - 2ai - 2D^ _ 2a2 + 2(1 - m)D^ + dA^ 

d + eA^ ~ b2-fA^ ■ ^'^^^ 

In the limit of m = 1, the periodic solution ()28p goes over to the hyperbolic dark-bright soliton solution 

(p = Atanh{Dx + xq) , 'ip = Bsech.{Dx + xq) , (38) 
satisfying the constraints ()35p to (j37p with m = 1. 



2.6 Solution VI 

Another dark-bright periodic soliton solution is given by 

(p = Ay/rnsn{Dx + xo,m) , = Bdn{Dx + xo,m) , (39) 



provided six coupled equations similar to (p9|) to (|Mj) are satisfied. Three of these equations determine the 
three unknowns A, B, D while the other three equations give three constraints between the nine parameters 
fli,2i &i,2j ci^2) d, e, f. For example, it is easily shown that two of the constraints are again given by Eq. ([36]) 
while and B^ are now given by 



A^ = ^ , B^ = ^ . (40) 

2ci 2c2 

while the other two constraints are 

2 _ biA^ - 2ai - 2mD'^ _ 2a2 - 2(1 - m)D^ + dA"^ 
~ d + eA^ ~ 62 - fA^ ■ ^ ^ 

In the limit of m = 1, the periodic solution (j39p goes over to the hyperbolic dark-bright soliton solution 
HMD. 



3 Solutions of Coupled 0^ Model In terms of Lame Polynomials of 
Order Two 

We now show that quite remarkably, the coupled model characterized by the field Eqs. ^ not only 
admits periodic solutions in terms of Lame polynomials of order one, but it also admits novel periodic 
solutions in terms of Lame polynomials of order two. It is worth reminding once again that neither Lame 
polynomials of order one nor of order two are solutions of the uncoupled (p^ problem. Since there are five 
Lame polynomials of order two, and since two of these are of the form ^sn^ [D{x + xq), m] -|- F, and further, 
the two field Eqs. ^ are symmetrical in (j) and ip, in principle there could be ten solutions of order two. 
However, it turns out that only two of these are admitted by the field Eqs. ^ which we now discuss. 
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3.1 Solution I 

It is easily shown that 

(p = Asn'^{Dx + xo,m) + F , ip = Bsn{Dx + xo,m)cn{Dx + xo,m) , (42) 

is an exact periodic solution to the coupled Eqs. ^ provided the following eleven coupled equations are 
satisfied 

aiF-biF^ + ciF^ = 2AD^ , (43) 

aiA - 36i^f2 + 5ciAF^ + eB'^F^ + dB^F = -4(1 + m)AD'^ , (44) 

- 3biA^F + lOciA^F^ + eB'^F'^{3A - F) + ^B^F + dB'^{A - F) = QmAD'^ , (45) 

f 

- biA^ + WciA^F^ + 3eAB^F{A - F) + -B^{A - 2F) - dAB"^ = , (46) 

bciA'^F + eA'^B'^{A - 3F) + ^B'^iF - 2A) = , (47) 

ciA'^ - eA'^B^ + ^B^ = . (48) 

a2 + ^F^ + dF^ = -{i + 771)0"^, (49) 

- b2B^ + 2eAF^ + fF^B^ + 2dAF = QmD'^ , (50) 

62 -B^ + C2S^ + SeA'^F^ + fB^F{2A - F) + dA^ = , (51) 

-2c2B'^ + 2eA^F + fAB'^{A-2F) =0, (52) 

C2B^ - fA^B^ + = . (53) 

Four of these equations determine the four unknowns A, B, D, F while the other equations give constraints 

between the nine parameters ai^2, ^1,2, ci^2,d, e, /. In particular, we find that the solution exists only if 

eA'' = fB\ e^ = 8clc2, f = 8clci , (54) 

and further if F 7^ 0. 

In the limit of m = 1, the periodic solution (|42p goes over to the hyperbolic solution 

(/) = Atanh^(Z)x + xo) + -F , -0 = tanh(L'x + 2;o)sech(-Dx + xq) , (55) 
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provided the constraints (03]) to (j53p with m = 1 are satisfied. There is one special case when this solution 
takes a simpler form, i.e. when A = —F, the solution is given by 

(j) = —Asech^{Dx + xq) , tfj = B tanh(Z)a; + XQ)sech{Dx + xq) , (56) 

provided Eq. ()54p is satisfied and further 

D'^ = ^ , ai = 4a2>0, d = -62 , bl{f - e) = 6a2C2e , 



B^ = '-^, ^.^±HhZ^. (57) 

62 ' 2ci ^ ^ 

Thus in the (f) variable, one is at T < since 6^ > Qa\C\. 

3.2 Solution II 

The other allowed solution is 

(j) = AsT?{Dx -\- XQ,m) + F ^ i/j = Bsn{Dx + xo,m)dn{Dx + xo,m) , (58) 

which is an exact periodic solution to the coupled Eqs. ([!]) provided Eqs. (H3]) . (fHI) and the following nine 
coupled equations are satisfied 

- SbiA'^F + lOciyl^F^ + eB'^F'^{3A - mF) + ^B^F + dB'^{A - mF) = 6mAD^ , (59) 

f 

- biA^ + WciA^F^ + 3eAB^F{A - mF) + ^B'^iA - 2mF) - dmAB^ = , (60) 

5ci^^F + eA'^B'^iA - 3mF) + ^5^(mF - 2A) = , (61) 

ci^^ - em^^S^ + ^m^S^ = . (62) 

ai + ^F^ + dF"^ = -{l + im)D\ (63) 

- 62^2 + 2eAF^ + /F^B^ + 2(i^F = 6mD'^ , (64) 
62m-B2 + 02^^ + 36^2^2 + fB^F{2A - mF) + dA^ = , (65) 

- 2mc25^ + 2eA^F + fAB'^{A - 2mF) = , (66) 

C2m^B^ - fmA'^B^ + = . (67) 
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Four of these equations determine the four unknowns A, B, D, F while the other equations give constraints 
between the nine parameters ai^2, ^1,2, ci^2,d, e, /. In particular, we find that the solution exists only if 

emA'^ = fB'^ , e' = 8cf C2 , f = 8c^ci , (68) 

and further if F 7^ 0. 

In the limit of m = 1, the periodic solution (ISSh also goes over to the hyperbolic solution (I55p . 

4 Coupled 06-04 Model 

We now consider a coupled 4>^-(l)^ model with bi-quadratic coupling. In particular, we consider the model 
characterized by the potential 

V{4>, ^) = 0^ - + + (1^^ + f + ^^^) + . (69) 

This leads to the coupled field equations 

Cpd) q o 

— = aic/) - hi(l)'^ + CKjf + #V , 

P^ = a2^ + /2^' + b2^' + d^<t? . (70) 

From stability considerations we shall always take ci > 0, 62 > 0. Note that in case /2 = 0, the model 
corresponds to the symmetric (/>^ model with a second order transition, while as long as /2 7^ 0, the model 
corresponds to a first order transition. We shall discuss the various solutions both when /2 7^ as well as 
when /2 = 0. 

4.1 Solution I 

It is easily shown that 



^Ia/I ± sn(L'x + xq, m) , ip = Bsti{Dx + XQ,m) + F , (71) 



is an exact periodic solution to the coupled Eqs. (|70|) provided the following seven coupled equations are 
satisfied 

n2 

ai-hiA^ + ciA^ + dF'^ = -— , (72) 
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- biA^ + 2ciA^ ± 2dBF = — ' (73) 

4 + dB^ = , (74) 

[a2 + f2F + b2F'^ + dA'^]F = , (75) 

+ 2f2FB + 362Si^^ + dA^{B ±F) = -{1 + m)BD^ , (76) 

[f2B + 3b2BF±dA'^]B = 0, (77) 

62-6^ = 2mL»2 . (78) 



Here, A and B denote the amplitudes of the kink lattice, D is an inverse characteristic length, F is a 
constant while xq is the (arbitrary) location of the kink. Four of these equations determine the four 
unknowns A, B, D, F while the other three equations give three constraints between the seven parameters 
ai,2,&i,2,ci,/2,d. 

In the limit of m = 1, the periodic solution ([7T]) goes over to the hyperbolic solution 

(j) = 1 ± tanh(i:>x + xq) , %l) = Biimh.{Dx + xq) + F , (79) 

provided Eqs. (1721) to ([78]) with m = 1 are satisfied. 
Several comments are in order at this stage. 

1. The solution (j7ip continues to exist even if /2 = 0, i.e. in the symmetric t/)^ case. 

2. Solution (I7ip continues to exist if F = 0. However, no solution exists in case i*" = /2 = as in that 
case d is also forced to be zero. 

3. In case B = itF then the solution exists only at m = 1. In particular it is easily shown that 

(/> = A^l±iimh.{Dx + XQ) , -0 = ^[1 ± tanh(i:>x + xq)] , (80) 

is an exact solution to field Eqs. (ffOj) provided 

D^ = a^, a2 = 4ai>0, ^2 = ^, B' = ^ . (81) 

01 ZO2 
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and further 

62ai(36f - 16aici) 62af (36f - 16aiCi) 

In the special case of /2 = 0, this solution continues to exist provided / = —3bi and hence (IGaiCi — 
3bj)aib2 = 6026?. 

4. On the other hand, the solution 



4> = A^l ± tanh(i:»j; + xq) , ip = F[l =F tanh(L>x + xq)] , (83) 
is an exact solution to field Eqs. (fTOl) provided 



D'^ + ai = biA^ , 2Z)2 + ai = 4ciA^ , - oi = 4^5^ , 

81)2 + 02 = 2/2^ , 41)2 - 02 = 2dA2 , 2L>2 = 62^^ • (84) 



These relations imply that AciA? = bi ± ^Jbf — 4aiCi. In the special case of /2 = 0, this solution 
exists provided a2 = —8D'^ , dA^ = 6D^ while the other four relations are as given by Eq. (j84p . 

4.2 Solution II 

We now present three solutions which are only valid at m = 1, i.e. in the hyperbolic limit. For example, 
it is easy to show that 

^ Asn{Dx + xo,m) Bcv?{Dx + XQ,m) 

~ 1 W — F TTF, i vT ' 



y/l- Fsn\Dx + XQ,m) ' [l-Fsn2(Da; + 

is an exact solution to the coupled Eqs. (|70p only if m = 1, and if 

ai + d-B^ = (3F - 2)Z)2 , bi = 2{l - F){D'^ + ai) , ciA^ = {I - FfilD"^ + ai) , 

a2 + f2B = -2{l + 3F)D\ dA^ = {1 - F)[6{1 + F)D^ - f2B] , b2B^ = 8D^F. (86) 

From here it follows that (bf — 4aiCi)>l'^ = 4(1 — F)'^D^ > 0. Note that at m = 1, the solution (fHSl) can be 
rewritten as 

^ ^tanh(Dx + xo) Bsech.'^{Dx + xq) 

' ^ 7i ttT 1 2/ r-i I vT ' V"' J 



^l-Ftanh2(i?x + xo) ' [I - F tanh^Dx + xq)] ' 

Note also that this solution continues to hold good even if /2 = 0. 
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4.3 Solution III 

Another solution, which is only valid at m = 1 is given by 



AsechiDx + xq) Bsech^iDx + xq) 



'l-Ftanh\Dx + xo) [1 " tanh {Dx + xq)] 

This is an exact solution to the coupled Eqs. ([70]) provided 

ai = D'^, 02 = 4ai , biA^ = 2(1 + F)D^ , 62^^ = 8FD^ , 
ciA"^ + dS^ ^ 3^2)2 , dyl^ + f^B = -6(1 + F)Z)2 . 

Note that this solution continues to hold good even if /2 = 0. 
4.4 Solution IV 

Yet another solution, which is only valid at m = 1 is given by 

A Bsech^Dx + xo) 

^ V = T-, — TTTT; — ] • (90) 



'1- Ftanh\Dx + xo) [1 " tanh {Dx + xq)] 

This is an exact solution to the coupled Eqs. ([70]) provided 

bi A"^ = 2{l - F){D^ + ai) , d-B^ + aiF^ = (3 - 2F)FD^ ciA"^ = {2D^ + ai){l - Ff , 

a2F + f2B = -2{3 + F)D^ , b2B^ = 8FD'^, dA^ = {1 - F){4D^ - a2) . (91) 

Prom here it follows that (6f — AaiCi)A'^ = 4(1 — F)'^D'^ > 0. Note that this solution continues to hold 
good even if /2 = 0. 

It is worth noting that for the three solutions as given by Eqs. ()87p . (I88p and ()90p . while (j) continues 
to be a solution of the uncoupled (j)^ field theory, in neither of these three cases, iIj is an exact solution of 
either the symmetric or the asymmetric, uncoupled ip^ problem. 

5 Conclusions 

In this paper we have shown that the Lame polynomials of order one and two are periodic solutions of a 
coupled (f)^ problem. These are novel solutions in the sense that while they are the solutions of the coupled 
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(f)^ problem, they are not the solutions of the corresponding uncoupled problems. In particular, we have 
obtained six solutions in terms of Lame polynomials of order one and two solutions in terms of Lame 
polynomials of order two. These results are applicable to both the structural phase transitions i2j and 
field theoretic contexts [3 SI [S] . 

We have also obtained four solutions of the coupled (f)^ — cf)^ problem, both when the c/)^ potential 
corresponds to a first order as well as a second order transition. Note that while the solutions of the coupled 
problem are also the solutions of the uncoupled (j)^ problem, but they are not the solutions of either the 
symmetric or the asymmetric uncoupled (/>^ problems. These solutions are also useful in understanding 
coexistence of different crystalline structures in elements \12\ [131 E! ^-^d ferroelectrics |151 [T6] . 

It will be interesting to obtain solutions of few other coupled field theories and with couplings that are 
not bi-quadratic. An example of a coupled model with linear-quadratic coupling occurs in the context of 
isostructural transitions [20j . It is conceivable that in some cases a linear-cubic coupling may be symmetry 
allowed. 
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